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ABSTRACT 

We  present  multifrequency  observations  with  the  Very  Long  Baseline  Array  of  the  compact  double 
radio  source  B2  2050  +  36.  Our  observations  are  at  0.33,  0.61,  1.67,  2.3,  and  8.4  GHz,  with  the  0.61  GHz 
observations  forming  the  third  epoch  of  observation  of  this  source  at  that  frequency.  At  0.61  GHz,  the 
structure  of  B2  2050  +  36  is  dominated  by  two  components  56  mas  apart.  Within  the  uncertainties  of  the 
various  measurements,  this  separation  has  remained  unchanged  for  the  past  16  yr.  Any  differential  image 
wander  caused  by  refractive  interstellar  scattering  is  less  than  4  mas.  Both  the  lack  of  differential  image 
wander  and  the  frequency  dependence  of  the  angular  diameter  of  B2  2050  +  36  below  1  GHz  indicate 
that  the  electron  density  power  spectrum  along  this  line  of  sight  has  a  spectral  index  near  the  Kolmogo¬ 
rov  value,  with  a  value  of  4  being  highly  unlikely.  We  conclude  that  diffractive  scattering  dominates 
along  this  line  of  sight;  results  in  the  literature  indicate  that  this  conclusion  also  holds  true  for  the  line  of 
sight  to  the  pulsar  PSR  B2020  +  28  (8?7  from  B2  2050  +  36).  Comparison  of  our  1.67  GHz  observations 
with  those  obtained  21  yr  previously  place  a  limit  on  the  projected  linear  separation  velocity  of  the  two 
components  of  c. 

Subject  headings :  galaxies :  active  —  galaxies :  individual  (B2  2050  +  36)  —  ISM :  structure  — 
radio  continuum :  ISM  —  scattering  —  turbulence 


1.  INTRODUCTION 

Density  fluctuations  on  AU-size  scales  in  the  interstellar 
plasma  may  produce  a  host  of  observable  effects  including 
refractive  intensity  scintillations  and  image  wander.  The 
magnitude  of  these  effects  depends  crucially  on  the  spec¬ 
trum  of  the  interstellar  density  fluctuations.  If  the  spectrum 
is  a  power  law  with  a  spectral  index  less  than  4,  these 
refractive  effects  will  generally  be  masked  by  diffractive 
effects  from  even  smaller  scale  density  fluctuations  (Romani, 
Narayan,  &  Blaifford  1986).  Conversely,  if  the  spectrum  is  a 
power  law  with  a  spectral  index  greater  than  4,  refractive 
effects  will  dominate.  For  many  lines  of  sight,  the  density 
spectral  index  appears  to  be  close  to,  but  less  than,  4 
(Armstrong,  Rickett,  &  Spangler  1995),  with  important 
exceptions  (see,  e.g.,  Hewish,  Wolszczan,  &  Graham  1985; 
Clegg,  Fiedler,  &  Cordes  1993;  Gupta,  Rickett,  &  Lyne 
1994;  Stinebring,  Faison,  &  McKinnon  1996;  Rickett, 
Lyne,  &  Gupta  1997). 

Complicating  the  effort  to  measure  refractive  effects  is 
their  generally  long  timescale.  The  relevant  timescale  is 
given  by  the  time  it  takes  for  an  AU-sized  density  fluctua¬ 
tion  to  drift  past  the  observer.  For  observations  of  pulsars, 
whose  velocities  can  exceed  100  km  s“^,  the  timescale  is 
days.  For  observations  of  extragalactic  sources,  the  relevant 
velocity  is  a  combination  of  the  Earth’s  velocity  and 
random  interstellar  motions  and  is  more  likely  to  be 
roughly  25  km  s“^,  meaning  that  the  refractive  timescale  is 
months  to  years. 

The  source  B2  2050  +  36  is  an  ideal  source  for  probing 
refractive  effects.  It  lies  on  the  outskirts  of  the  Cygnus 
superbubble  (Bochkarev  &  Sitnik  1985),  a  region  known  to 
exhibit  enhanced  interstellar  scattering  (Fey,  Spangler,  & 
Mutel  1989;  Lazio,  Spangler,  &  Cordes  1990;  Fey,  Span¬ 


gler,  &  Cordes  1991;  Wilkinson,  Narayan,  &  Spencer  1994; 
Molnar  et  al.  1995;  Spangler  &  Cordes  1998).  VLBI  obser¬ 
vations  reveal  its  milliarcsecond  structure  to  be  dominated 
by  two  components  separated  by  approximately  60  mas; 
thus,  a  search  for  refractive  image  wander  requires  only 
relative  position  measurements. 

Fey  &  Mutel  (1993,  hereafter  FM93)  obtained  0.61  GHz 
VLBI  observations  of  B2  2050  +  36.  Combined  with  pre¬ 
vious  measurements  in  the  literature,  they  showed  that  scat¬ 
tering  is  moderately  strong  along  this  line  of  sight.  They 
were  also  able  to  constrain  the  amount  of  differential 
refractive  image  wander  between  the  two  components. 
Their  constraint,  however,  is  based  on  measurements  of  the 
separation  of  the  components  at  two  epochs  separated  by  4 
yr.  If  the  refractive  timescale  is  longer  than  4  yr,  their  con¬ 
straint  is  weakened  considerably.  Assuming  that  typical 
interstellar  velocities  along  this  line  of  sight  are  50-100  km 
s“^,  FM93  argued  that  the  refractive  timescale  is  compara¬ 
ble  to  or  less  than  4  yr.  Obviously,  if  the  typical  interstellar 
velocity  is  more  like  25  km  s“^,  the  refractive  timescale  is 
much  longer. 

We  have  acquired  a  third  epoch  of  0.61  GHz  VLBI  obser¬ 
vations  of  B2  2050  +  36,  specifically  to  confront  the  conclu¬ 
sions  of  FM93  with  another  epoch  of  observations.  This 
third  epoch  is  12-16  yr  after  the  first  two  epochs  and  should 
allow  us  to  place  much  better  constraints  on  any  refractive 
wander.  We  have  also  acquired  simultaneous  or  quasi- 
simultaneous  observations  at  a  range  of  other  wavelengths 
in  order  to  confirm  previous  measurements  of  the  strength 
of  interstellar  scattering  along  this  line  of  sight. 

In  §  2  we  describe  our  observations  and  data  analysis,  in 
§  3  we  discuss  the  implications  of  our  observations  for  both 
diffractive  and  refractive  interstellar  scattering  along  this 
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TABLE  1 

VLBA  Observing  Log 


Frequency  Bandwidth  Recorded  On-source  Time  Synthesized  Beam  Image  Noise  Level* 
(GHz)  (MHz)  Polarization  (hr)  (mas)  (mjybeam"*) 


0.33 .  12  R,  L  2.88”  61  x  58  at  87°  3.2 

0.61  .  4  R,  L  2.88”  27  x  19  at  14°  5 

1.67  .  32  R,  L  2.93  7.3  x  4.7  at  -T.9  0.78 

2.27 .  16  R,  L  2.92°  5.8  x  3.6  at  -9?4  1.1 

8.42 .  16  R,  L  2.92°  1.9  x  1.4  at -10°  0.61 


*  All  images  are  of  Stokes  I  polarization. 

”  Observations  at  0.33  and  0.61  GHz  were  conducted  simultaneously. 
°  Observations  at  2.27  and  8.42  GHz  were  conducted  simultaneously. 


line  of  sight,  and  in  §  4  we  discuss  the  limits  our  obser¬ 
vations  place  on  the  intrinsic  structure  and  kinematics  of 
the  source.  We  present  our  conclusions  in  §  5. 

2.  OBSERVATIONS  AND  DATA  ANALYSIS 

We  observed  B2  2050  +  36  with  the  Very  Long  Baseline 
Array  (VLBA)  for  a  total  of  10  hr  on  2000  July  25.  Relevant 
observation  parameters  are  summarized  in  Table  1. 

The  observation  frequencies  were  0.33, 0.61, 1.67,  2.3,  and 
8.4  GHz;  both  right  and  left  circular  polarizations  were 
recorded  at  all  frequencies.  The  2.3  and  8.4  GHz  obser¬ 
vations  were  obtained  simultaneously  in  “S/X  mode.” 
Simultaneous  0.33  and  0.61  GHz  observations  also  were 
obtained  by  setting  two  of  the  eight  baseband  converters 
(BBCs)  at  each  VLBA  antenna  to  record  the  two  senses  of 
circular  polarization  at  0.61  GHz.  Because  of  severe  radio 

TABLE  2 


Gaussian  Source  Models 


V 

(GHz) 

S 

(Jy) 

r 

(mas) 

0 

(deg) 

a 

(mas) 

bj  a 

<!> 

(deg)  Component 

0.33.... 

2.17 

0.0 

106.2 

0.76 

-35 

E* 

0.61.... 

2.48 

0.0 

24.1 

0.82 

24 

E 

0.39 

56.1 

-82 

26.2 

0.69 

8 

W 

0.11 

32.6 

8 

47.5 

0.00 

67 

E 

1.67.... 

1.75 

59.3 

-80 

3.5 

0.90 

-49 

W 

0.79 

5.5 

-82 

4.2 

0.74 

-19 

E 

0.99 

1.3 

-132 

4.0 

0.86 

33 

E 

1.06 

3.4 

42 

4.9 

0.69 

-77 

E 

0.16 

8.9 

-106 

11.8 

0.32 

-52 

E 

2.27.... 

1.97 

0.0 

-32 

2.9 

0.71 

-61 

W 

0.92 

58.5 

101 

3.7 

0.81 

17 

E 

0.58 

53.8 

100 

3.3 

0.59 

-25 

E 

0.83 

61.5 

97 

4.2 

0.54 

-78 

E 

0.12 

52.9 

103 

7.5 

0.39 

-46 

E 

8.42.... 

0.35 

18.4 

-89 

0.8 

0.54 

84 

W 

0.39 

20.1 

-87 

1.0 

0.51 

-1 

W 

0.19 

34.9 

103 

2.1 

0.64 

-29 

E 

0.19 

40.3 

105 

2.1 

0.54 

60 

E 

0.19 

43.2 

99 

3.5 

0.51 

-85 

E 

0.22 

18.8 

-85 

3.1 

0.24 

-67 

W 

0.16 

39.0 

102.4 

6.3 

0.28 

40 

E 

Note. — The  models  fitted  to  the  visibility  data  are  Gaussians  with  flux 
density  S  and  FWHM  major  axis  a  and  minor  axis  6,  with  major  axis  in 
position  angle  (f)  (measured  north  through  east).  Components  are  separat¬ 
ed  from  the  (arbitrary)  origin  of  the  image  by  a  distance  r  in  position  angle 
6,  which  is  the  position  angle  (measured  north  through  east)  of  a  line 
joining  the  components  with  the  origin.  Components  and  subcomponents 
are  identified  as  with  either  an  “E”  or  “W”  corresponding  to  whether 
they  are  eastern  or  western,  respectively. 

*  Assumed  eastern  component;  see  §  3.1. 


frequency  interference  around  0.61  GHz,  filters  restricted 
the  observable  bandwidth  for  each  BBC  to  4  MHz.  The 
remaining  BBCs  recorded  at  0.33  GHz.  Formatter  prob¬ 
lems  at  the  Kitt  Peak  antenna  meant  that  only  nine  of  the 
10  VLBA  antennas  were  correlated. 

The  correlated  visibilities  were  processed  in  the  standard 
fashion  within  the  NRAO  AIPS.  Fringe-fitting  intervals 
ranged  from  10  minutes  at  1.67  GHz  to  3  minutes  at  0.33 
and  8.4  GHz.  At  0.33  and  0.61  GHz,  substantial  editing  of 
the  calibrated  visibilities  was  required  to  remove  radio  fre¬ 
quency  interference  (RFI). 

In  general,  the  calibrated  visibilities  were  then  exported 
to  the  Caltech  DIFMAP  program  for  hybrid  mapping.  To 
ensure  a  robust  estimate  for  the  separation  of  the  com¬ 
ponents  (§  3.2),  the  0.61  GHz  observations  were  imaged 
using  both  AIPS  and  DIFMAP;  the  0.33  GHz  observations 
were  imaged  using  only  AIPS. 

Figure  1  shows  the  images  of  B2  2050  +  36  at  the  various 
frequencies.  The  observed  rms  noise  level  in  these  images 
ranges  from  0.61  mJy  beam”^  at  8.42  GHz  to  5  mJy 
beam”  ^  at  0.61  GHz  and  is  summarized  in  Table  1.  The  rms 
noise  levels  in  our  images  are  roughly  within  a  factor  of  2-3 
of  the  expected  value. 

Using  the  images  shown  in  Figure  1  as  a  guide,  we  fitted 
the  self-calibrated  visibility  data  with  multiple  Gaussian 
component  models.  The  number  of  Gaussian  components 
fitted  ranged  from  1  at  0.33  GHz  to  7  at  8.42  GHz.  These 
models  are  summarized  in  Table  2.  Below  we  shall  make  use 
of  the  fitted  angular  diameters.  At  the  higher  frequencies,  we 
take  the  uncertainty  in  a  fitted  angular  diameter  to  be  one- 
half  of  the  synthesized  beam;  at  the  lower  frequencies, 
where  a  one-  or  two-component  fit  suffices,  we  take  the 
uncertainties  to  be  the  square  roots  of  the  diagonal  elements 
of  the  covariance  matrix  obtained  from  a  least-squares  fit  to 
the  visibility  data  (provided  by  the  UVFIT  fitting  routine 
within  AIPS);  these  values  are  on  the  order  of  0.5  mas. 

3.  CONSTRAINTS  ON  INTERSTELLAR  SCATTERING 

We  now  revisit  the  scattering  properties  of  the  line  of 
sight  to  B2  2050  +  36.  We  consider  both  diffractive  effects 
(namely,  angular  broadening)  and  refractive  effects  (namely, 
differential  image  wander  and  refractive  intensity 
fluctuations). 

We  shall  make  the  standard  assumption  throughout  that 
the  density  fluctuations  responsible  for  interstellar  scat¬ 
tering  can  be  parameterized  by  an  exponentially  truncated 
power  law 

=  q>  q^  .  (1) 
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Fig.  1. — VLSI  images  of  B2  2050  +  36.  (a)  0.33  GHz.  The  rms  noise  level  is  3.2  mjy  beam"  \  and  the  contour  levels  are  (3.2  x  —2,  3,  5,  7.1,  10 ... )  mjy 
beam"\  The  beam  is  60.8  x  58.5  mas  at  a  position  angle  of  87°  and  is  shown  in  the  lower  left-hand  corner,  (b)  0.61  GHz.  The  rms  noise  level  is  5  mJy  beam"*, 
and  the  contour  levels  are  (5  x  —3,  3,  5,  7.1,  10...)  mJy  beam"*.  The  beam  is  27  x  19  mas  at  a  position  angle  of  14°  and  is  shown  in  the  lower  left-hand 
corner,  (c)  1.67  GHz.  The  rms  noise  level  is  0.78  mjy  beam"*,  and  the  contour  levels  are  (0.78  x  —3,  3,  5, 7.1, 10. . .)  mJy  beam"*.  The  beam  is  7.3  x  4.7  mas 
at  a  position  angle  of  —  7?9  and  is  shown  in  the  lower  left-hand  corner,  (d)  2.3  GHz.  The  rms  noise  level  is  1.1  mJy  beam"*,  and  the  contour  levels  are 
(1.1  X  —3,  3,  5,  7.1, 10. ..)  mJy  beam"*.  The  beam  is  5.8  x  3.6  mas  at  a  position  angle  of  —  9?4  and  is  shown  in  the  lower  left-hand  corner,  (e)  8.4  GHz.  The 
rms  noise  level  is  0.61  mJy  beam"*,  and  the  contour  levels  are  (0.61  x  —3,  3,  5,  7.1,  10...)  mJy  beam"*.  The  beam  is  1.9  x  1.4  mas  at  a  position  angle  of 
—  10°  and  is  shown  in  the  lower  left-hand  corner. 


Here  is  the  largest  wavenumber  on  which  density  fluctua¬ 
tions  occur,  which  is  related  to  the  smallest  length  scales  for 
density  fluctuations  or  the  “inner  scale,”  as 
The  smallest  wavenumber  on  which  density  fluctuations 
occur,  or  the  “  outer  scale,”  is  q^  =  InflQ.  The  constant  is 
a  measure  of  the  level  of  density  fluctuations  and  is  related 


to  the  rms  density  as  {nj}  =  I2(2n)'^~°‘/(a  — 

(Cordes  et  al.  1991).  The  power-law  spectral  index  is 
assumed  to  provide  clues  about  the  process(es)  that  gener¬ 
ate  or  maintain  the  density  fluctuations.  Our  observations 
allow  us  to  measure  a,  constrain  C^,  and  indirectly  place 
limits  on  . 
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Fig.  le 


3.1.  Diffractive  Scattering 
3.1.1.  Angular  Broadening  and  Density  Spectrum 

An  infinitely  distant  point  source  viewed  through  a 
medium  filled  with  density  fluctuations  having  a  spectrum 
described  by  equation  (1)  and  a  spectral  index  a  =  11/3  will 
be  broadened  to  an  angular  diameter  (Cordes  &  Lazio 
1991 ;  Taylor  &  Cordes  1993) 

'128  mas  , 

=  J  /I 

1l32masSMi/SGi(^^^^j  ,  0,  >  . 

(2) 

Here  Vqhz  is  the  observing  frequency  measured  in  GHz  and 
the  scattering  measure  SM  is  the  line  of  sight  integral  of  C^, 

SM  =  I  dz  C„^(z)  .  (3) 

The  quantity  is  the  crossover  angle  and  describes  the 
relative  importance  of  the  inner  scale  It  is  given  by 
Across  —  0"16vdHz(ii/100  km)“^.  The  constants  (128  and  132 
mas)  given  in  equation  (2)  are  those  appropriate  for  infi¬ 
nitely  distant  sources  and  differ  slightly  from  those  given  by 
Cordes  &  Lazio  (1991),  who  considered  sources  embedded 
in  the  medium.  These  formulae  specify  the  diameter  as 
FWHM  of  the  shape  of  the  image,  and  that  will  be  the 
convention  we  adopt  in  reporting  the  measured  diameters. 

The  scattering  diameter  6^  defines  the  diffractive  scale 
frf  =  i/drf  (at  observing  wavelength  X).  If  the  typical  interfer¬ 
ometer  baseline  b  is  sufficient  to  resolve  6^  (i.e.,  b  ^  l^)  and 
b  <li,  we  expect  6^  >  and  0^oc  v“^.  Conversely  if 
b  >  fi,  then  we  expect  6^  <  6^^^^  and  6^  oc  There¬ 

fore,  fitting  the  observed  diameter  as  a  function  of  frequency 
can  constrain  both  and  a. 

There  are  alternate  causes  for  a  6^az\~^  dependence. 
For  instance,  if  the  spectrum  of  density  fluctuations  is  not  a 
power  law  but  a  Gaussian  centered  on  a  particular  length 


scale,  that  can  produce  a  0^  oc  dependence.  Similarly,  if 
the  density  spectrum  is  dominated  just  slightly  by  AU-size 
scales  such  that  a  >  4,  that  could  also  produce  an  apparent 
0^  oc  dependence  within  the  measurement  uncertainties. 
As  we  shall  show  below,  a  <  4  so  that  we  could  interpret  a 
0^  oc  dependence  as  resulting  from  a  “large”  inner 
scale. 

Figure  2  shows  the  measured  angular  diameters  as  a 
function  of  observing  frequency.  Also  shown  are  two  curves, 
one  for  0  oc  V  appropriate  for  a  spectrum  with  a  =  11/3 
and  a  diffractive  scale  much  larger  than  and  one  for 
0  oc  appropriate  for  the  case  in  which  h^h- 

As  Figure  1  demonstrates,  B2  2050+36  is  not  only  not  a 
point  source,  but  the  two  components  of  it  can  be  resolved 
into  multiple  subcomponents  at  high  frequencies.  A  visual 
inspection  of  Figure  2  suggests  that  at  frequencies  below  1 
GHz,  scattering  dominates  the  structure  of  B2  2050  +  36, 
while  at  frequencies  above  1  GHz,  there  are  clear  deviations 
from  either  curve  presumably  due  to  intrinsic  structure.  In 
order  to  quantify  this  result,  we  have  fitted  the  measured 
diameters  of  the  eastern  component  (Table  2)  with  the  fol¬ 
lowing  functional  form : 

0^  =  0?  +  .  (4) 

Here  0;  is  the  intrinsic  diameter  and  0^  is  the  scattering 
diameter,  both  at  the  fiducial  frequency  of  1  GHz.  Unlike  in 
equation  (2),  in  which  we  specified  the  dependence  for  the 
specific  case  of  a  =  1 1/3,  here  we  allow  the  frequency  depen¬ 
dence  to  be  part  of  the  lit.  The  frequency  dependence  we  use 
is  appropriate  for  a  <  4. 

Several  points  are  in  order  regarding  the  choice  of  data 
we  used.  At  0.33  GHz,  only  one  component  is  seen  (Fig.  la). 
Attempts  to  fit  the  visibility  data  with  a  two-component 
structure  suggested  by  the  0.61  GHz  image  produced  a 
poorer  fit  to  the  data.  We  have  subtracted  a  one-component 
model  (Table  2)  from  the  visibility  data  and  imaged  the 
residuals.  There  is  no  indication  of  a  second  component  at 
the  location  expected  from  the  higher  frequency  images.  We 
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Fig.  2. — Angular  diameter  of  the  components  of  B2  2050  +  36  as  a  function  of  observing  frequency.  All  components  listed  in  Table  2  are  plotted.  Top 
panel:  The  eastern  component.  Bottom  panel:  The  western  component.  For  the  western  component,  we  also  show,  as  an  open  circle,  the  diameter  it  would 
have  assuming  if  its  scattering  diameter  were  identical  to  that  of  the  eastern  component.  This  point  is  shown  for  display  purposes  only  and  is  not  used  in  any 
of  the  fits  described.  In  both  panels  a  dashed  line  shows  the  dependence  0  oc  v “ ^  as  expected  for  a  Kolmogorov  spectrum,  and  a  dotted  line  shows  8  ocv~^. 
Both  lines  are  normalized  to  pass  through  the  measured  angular  diameter  at  0.61  GHz.  For  the  eastern  component  only,  a  solid  line  shows  the  fit,  which 
allows  for  intrinsic  structure  and  which  is  discussed  in  the  text.  At  the  higher  frequencies,  we  plot  all  subcomponents,  but  we  show  the  uncertainty  on  only  one 
of  them  since  we  assume  the  uncertainty  to  be  the  same  for  all  subcomponents  and  equal  to  one-half  of  the  synthesized  beam. 


have  also  produced  model  visibility  data,  comprised  of  two 
components  whose  flux  densities,  diameters,  and  positions 
are  consistent  with  those  seen  at  higher  frequencies.  Fitting 
these  model  data  with  a  one-component  model  results  in  a 
fit  nearly  the  same  as  what  we  And  by  fitting  the  observed 
visibility  data.  The  similarity  suggests  that  any  systematic 
error  in  determining  the  diameter  is  probably  no  more  than 
1  mas.  Because  of  the  self-calibration  we  performed,  we 
have  no  absolute  position  information.  Given  that  the 
eastern  component  is  stronger  at  0.61  GHz,  we  have  there¬ 
fore  assumed  that  the  single  component  seen  at  0.33  GHz  is 
the  eastern  component. 

At  the  higher  frequencies,  both  components  can  be  re¬ 
solved  into  multiple  subcomponents.  In  Figure  2  we  plot 
the  diameters  of  the  individual  subcomponents.  The  diam¬ 
eters  of  the  individual  subcomponents  are  clearly  not  repre¬ 
sentative  of  the  diameter  of  the  entire  component. 
Therefore,  as  a  measure  of  the  size  of  the  components  at 
these  frequencies,  we  have  added  the  diameters  of  the  sub¬ 
components  (Table  2)  in  quadrature.  If  anything,  this  esti¬ 
mate  will  underestimate  the  size  of  a  component.  However, 
our  primary  interest  is  not  in  the  intrinsic  structure,  and  this 
procedure  suffices  to  obtain  an  estimate  for  the  intrinsic 
size. 

B2  2050  -I-  36  is  classified  as  a  gigahertz-peaked  spectrum 
source,  and  both  components  show  a  single  flux  density 
maximum  near  1  GHz  (Mutel,  Hodges,  &  Phillips  1985). 
Accordingly,  in  equation  (4)  we  have  assumed  that  the 
intrinsic  structure  has  no  frequency  dependence.  While  we 


have  not  conducted  extensive  tests,  assuming  no  frequency 
dependence  produces  clearly  much  better  fits  than  does 
assuming  a  frequency  dependence  of  v“^,  as  might  be 
expected  for  a  homogeneous  source  with  a  peak  brightness 
temperature  7^  (Kellermann  &  Owen  1988). 

We  have  used  a  grid-search  technique  to  minimize  the 
difference,  in  a  sense,  between  the  data  and  the  expres¬ 
sion  of  equation  (4)  for  the  three  parameters :  0;,  6^,  and  a. 
We  have  also  determined  the  allowed  ranges  for  the  param¬ 
eters  in  the  following  manner:  holding  two  parameters  fixed 
at  the  values  that  minimize  x^,  we  change  the  third  param¬ 
eter  until  x^  changes  by  unity. 

Table  3  summarizes  the  best-fit  value  for  the  parameters 
and  their  allowed  ranges.  The  best-fitting  curve  is  also 
shown  in  Figure  2.  A  key  result  of  this  fit  is  that  a  density 


TABLE  3 

Best-Fit  Diffractive  Scattering  Model  Parameters 
FOR  Eastern  Component  at  1  GHz 


Parameter 

Adopted  Value 

Range 

0j  (mas) . 

3.9 

3.5-4.3 

0,  (mas) . 

7.0 

6.98-7.03 

a . 

3.52 

3.51-3.53 

Derived  Scattering  Parameters 

SM‘‘(kpcm-^'>'5) . 

io-^‘ 

li  (km) . 

■^2000 

“  Kolmogorov  density  spectrum,  a  =  11/3,  assumed. 
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spectrum  spectral  index  of  4  is  highly  unlikely.  This  result 
means  that  the  inner  scale  is  substantially  smaller  than 
most  baselines  in  the  array  or  that  6^  <  ^^ross-  The  typical 
array  baseline  is  roughly  2000  km,  so  we  place  an  upper 
limit  of  1 1  2000  km.  This  constraint  is  not  a  serious  chal¬ 

lenge  to  the  value  estimated  by  Spangler  &  Gwinn  (1990)  of 
li  <  200  km.  Using  equation  (2),  we  can  also  derive  the 
scattering  measure  SM.  Table  3  also  lists  the  values  of  these 
derived  quantities. 

3.1.2.  Image  Anisotropy  and  Orientation 

The  anisotropy  and  orientation  of  a  scattered  image  can, 
in  principle,  yield  information  about  the  density  structures 
responsible  for  the  angular  broadening  (see,  e.g.,  Wilkinson 
et  al.  1994;  Trotter,  Moran,  &  Rodriguez  1998;  Cordes  & 
Lazio  2001).  The  simultaneous  low-frequency  observations 
we  report  here  are  seemingly  ideal  for  constraining  any 
wavelength-dependent  changes  in  the  anisotropy  or  orien¬ 
tation  of  B2  2050  +  36  since  the  u-v  plane  coverage  should 
be  identical  (once  scaled  by  the  ratio  of  the  wavelengths). 

We  consider  first  the  anisotropy  of  the  image.  We  find 
axial  ratios  of  b/a  ^  0.80  (Table  2)  or,  in  terms  of  the  image 
anisotropy  parameter  of  Romani  et  al.  (1986),  ^  0.25. 

Romani  et  al.  (1986)  also  provide  an  estimate  of  the  image 
anisotropy  to  be  expected  from  random  fluctuations  within 
the  medium.  For  our  observing  wavelengths  and  scattering 
parameters  derived  above,  the  expected  random  contribu¬ 
tion  to  is  approximately  0.03.  We  conclude  that  these 
image  anisotropies  are  significant.  These  anisotropies  are 
also  comparable  to  those  of  other  sources  seen  through  the 
Cygnus  region  (Spangler  &  Cordes  1988,  1998;  Wilkinson 
et  al.  1994;  Molnar  et  al.  1995;  Desai  &  Fey  2001).  Desai  & 
Fey  (2001)  attribute  the  similarity  of  the  axial  ratios 
between  different  sources  at  different  frequencies  and  differ¬ 
ent  epochs  to  be  indicative  of  an  anisotropy  in  the  scat¬ 
tering  material  itself. 

We  cannot  conclude  that  the  image  anisotropy  changes 
with  frequency  for  B2  2050+36.  The  formal  uncertainties 
from  the  fitting  of  the  visibility  data  are  approximately  1%. 
Even  this  small  uncertainty  on  the  angular  diameters,  when 
propagated  to  the  axial  ratios,  implies  an  uncertainty  of 
A(b/a)  ^  0.03,  comparable  to  the  difference  between  the 
axial  ratios  at  the  two  frequencies  (Table  2). 

Changes  in  the  image  orientation  between  the  two  fre¬ 
quencies  may  reflect  changes  in  the  orientation  of  the 
density  irregularities  responsible  for  the  scattering.  Unfor¬ 
tunately,  we  also  cannot  draw  any  conclusions  about  the 
image  orientation.  Although  the  u-v  plane  coverage  should 
be  similar,  differences  in  the  editing  (due  to  the  significant 
levels  of  RFI  encountered)  change  the  relative  u-v  plane 
coverage  significantly.  An  indication  of  the  difference  can  be 
seen  by  inspecting  the  synthesized  beams  shown  in  Figures 
la  and  lb.  In  particular,  the  synthesized  beam  at  610  MHz 
has  an  orientation  of  14°  compared  to  an  orientation  for  the 
eastern  component  of  24°.  Thus,  the  apparent  change  in 
orientation  between  330  (—35°)  and  610  MHz  (24°)  may  be 
largely  due  to  the  differences  in  u-v  coverage. 

3.2.  Refractive  Interstellar  Scattering 

Large-scale  density  fluctuations  may  affect  the  image 

position  and  flux  density.  In  this  section  we  place  con¬ 
straints  on  the  strength  of  refractive  modulation  of  the 
image,  focusing  primarily  on  the  0.61  MHz  image. 


3.2.1.  Differential  Image  Wander 

Density  fluctuations  on  a  scale  l^  =  D6^,  where  D  is  the 
characteristic  distance  to  the  scattering  medium,  will  cause 
refractive  “bending”  and  “steering”  of  the  propagating 
rays.  In  turn,  this  will  affect  the  angle  of  arrival  of  rays  and 
make  radio  sources  wander  about  a  true  or  nominal  posi¬ 
tion.  The  magnitude  of  this  effect  depends  both  on  the 
strength  of  scattering  and  the  value  of  a. 

FM93  selected  B2  2050  +  36  because  its  two-component 
structure  allowed  for  relative  position  measurements  to  be 
performed,  obviating  the  need  for  determining  absolute 
positions  with  milliarcsecond  accuracy.  FM93  summarized 
the  predictions  for  the  amount  of  differential  angular 
wander  as  a  function  of  a  as  predicted  by  various  groups 
(Cordes,  Pidwerbedsky,  &  Lovelace  1986;  Romani  et  al. 
1986;  Rickett  &  Coles  1988;  see  Appendix  A  of  FM93). 
They  found  that  the  rms  differential  refractive  angular 
wander  is 


(a -4)1 2 

«<4, 

-  ■'  ly  \(0i-4)l2  (5) 

G(a)(^^  j  0,  ,  a  >  4  . 

Here  is  the  separation  between  pierce  points  through 
the  scattering  medium.  For  a  double  source  (like  B2 
2050  +  36)  with  a  component  angular  separation  of  ^^^p, 
Gep  =  ^^sep-  The  functions  F{a)  and  G(a)  are  normalization 
functions  that  differ  from  group  to  group  but  typically  by 
no  more  than  a  factor  of  a  few.  For  reference,  Cordes  et  al. 
(1986)  predict  =  0.18  mas{Dy.pJ2.^^~ for  a 

Kolmogorov  spectrum  in  a  screen  located  fj^pp  kpc  distant, 
producing  an  angular  broadening  of  6^  mas,  and  observed 
at  a  wavelength  of  2.^^  cm. 

In  general,  if  a  >  4,  refractive  effects  like  angular  wander 
dominate  over  diffractive  effects  like  angular  broadening. 
Thus,  one  would  expect  to  see  the  apparent  position  of  a 
source  wander  by  more  than  the  width  of  its  scattering  disk. 
If  a  <  4,  the  converse  is  true. 

Table  4  summarizes  the  0.61  MHz  measurements  of  the 
angular  separation  between  the  two  components  of  B2 
2050  +  36.  We  have  estimated  using  two  different 

methods.  Henceforth,  we  shall  also  distinguish  between  the 
ensemble  average  expected  differential  angular  wander 
<A0r  and  the  quantity  estimated  from  the  separation 
measurements  First,  we  calculated  using  the 

standard  estimator  for  the  variance  and  using  the  three 
values  in  Table  4.  This  estimate  is  (0  ^  )^^^  =  2.3  mas.  Second, 
we  estimated  to  be  comparable  to  the  maximum 


TABLE  4 

B2  2050  +  36  Component  Separation  Measurements 


Epoch  (mas)  Reference 

1984  Aug  .  60  ±  3  Mutel  &  Hodges  1986 

1988  Jun .  56  ±  2  FM93 

2000  Jul .  56.1  ±  1.0  This  work 

Differential  Angular  Wander  Upper  Limits 

(mas) .  2.3 

max  (A0)  (mas) .  4 
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difference  between  component  separations  at  the  different 
epochs  max  (A6).  This  estimate  is  max  (A0)  =  4  mas. 

Figure  3  shows  as  a  function  of  the  density 

spectral  index  a,  as  predicted  by  Cordes  et  al.  (1986), 
Romani  et  al.  (1986),  and  Rickett  &  Coles  (1988).  We  also 
plot  our  observational  limits  {6  3(0  and  max  (A0). 

In  producing  Figure  3,  we  have  assumed  D  =  2  kpc.  We 
have  considered  <A0^>^^^  as  a  function  of  both  D  and  a. 
However,  the  dependence  on  D  is  fairly  weak,  oc 

for  a  =  11/3.  Reasonable  ranges  for  D  produce  little 
difference  in  The  value  D  =  2  kpc  is  motivated  by 

the  distance  to  the  Cygnus  superbubble  (Bochkarev  & 
Sitnik  1985),  through  which  the  line  of  sight  to  B2  2050  +  36 
passes. 

Figure  3  shows  clearly  that  along  this  line  of  sight  the 
upper  limits  on  constrain  a  <  4.  The  largest 

allowed  value,  a  =  3.86,  occurs  if  we  adopt  an  upper  limit  of 
3(0  jY^^  and  use  the  expression  for  <A0/>^^^  from  Rickett  & 
Coles  (1988).  Using  a  more  stringent  limit  on  the  differential 
angular  wander  max  (A0),  the  expressions  for  <A0/>^^^ 
from  Cordes  et  al.  (1986)  or  Romani  et  al.  (1986),  we  find 
a  <  3.7.  We  adopt  a  =  11/3  as  a  limit  based  on  (0/)^^^  we 
consider  to  be  overly  strict. 


3.6  3.8  4  4.2  4.4 

POWER  LAW  SPECTRAL  INDEX  a 


Fig.  3. — Differential  angular  wander  as  a  function  of  the 

density  spectral  index  a.  Three  expressions  for  are  shown:  that  of 

Cordes  et  al.  (1986;  dot-dashed  line),  Romani  et  al.  (1986;  dotted  line),  and 
Rickett  &  Coles  (1988;  triple-dot-dashed  line).  All  expressions  are  taken 
from  Appendix  A  of  FM93.  The  abrupt  change  in  slope  at  a  =  4  oceurs 
because  no  attempt  has  been  made  to  join  smoothly  the  expressions  for 
a  <  4  and  a  >  4.  The  solid  horizontal  lines  represent  our  upper  limits  on 
the  differential  angular  wander  of  B2  2050  +  36  as  given  by  the  standard 
variance  estimator  of  (6^)^'^  and  3(6 The  dashed  horizontal  line  rep¬ 
resents  our  upper  limit  on  the  differential  angular  wander  as  given  by 
max  (A0). 


The  standard  estimator  for  the  variance,  which  we  have 
used  in  determining  (0  ^ )^^^,  assumes  that  the  measurements 
are  independent.  In  this  case,  the  independence  of  the 
various  measurements  can  be  assessed  by  comparing  the 
interval  between  the  different  measurement  epochs  with  the 
refractive  timescale. 


20  yr 


D 

2  kpc 


10  km  s“ 


(6) 


where  v  is  the  bulk  velocity  of  the  material.  FM93  assumed 
V  =  50-100  km  s“^,  for  which  At^  =  2-4  yr.  In  this  case,  all 
three  measurement  epochs  (Table  4)  could  be  considered 
independent,  although  the  first  two  only  marginally  so. 
However,  if  the  characteristic  velocity  of  the  scattering 
medium  is  smaller,  by  as  httle  as  a  factor  of  2,  the  refractive 
timescale  would  be  long  enough  that  the  first  two  epochs 
(1984  August  and  1988  June)  could  not  be  considered  inde¬ 
pendent. 

If  the  characteristic  velocity  is  as  small  as  10  km  s“  then 
less  than  one  refractive  time  has  elapsed  since  the  first 
observation,  and  none  of  these  observations  can  be  con¬ 
sidered  to  be  independent.  Even  if  this  was  the  case,  we 
would  expect  to  see  some  change  over  the  12  yr  time  span 
separating  the  second  and  third  epochs  (1988  June  and  2000 
July).  That  essentially  no  change  in  the  angular  separation 
is  seen  indicates  that  max  (A0)  provides  a  reasonably  strin¬ 
gent  upper  limit  on  <A0/>^^^. 


3.2.2.  Refractive  Intensity  Fluctuations 

In  concert  with  refractive  wandering,  the  large-scale 
density  fiuctuations  will  also  produce  intensity  fiuctuations 
as  more  or  fewer  rays  are  “  steered  ”  into  the  line  of  sight.  In 
the  case  of  B2  2050  +  36,  this  will  produce  independent 
variations  in  the  fiux  density  of  its  two  components. 

FM93  compared  the  fiux  densities  of  the  two  components 
between  their  observations  and  those  of  Mutel  &  Hodges 
(1986).  They  found  that  the  flux  density  of  the  components 
was  consistent  with  no  change  in  the  4  yr  between  1984 
August  and  1988  June.  Table  5  summarizes  these  earlier 
measurements  and  those  reported  here. 

Because  of  changes  in  telescopes,  receivers,  correlators, 
and  possible  intrinsic  changes,  we  do  not  compare  the  abso¬ 
lute  flux  densities  of  these  earlier  epochs  to  our  measure¬ 
ments.  Rather,  we  compare  the  flux  density  ratio  between 
the  two  components.  As  Table  5  shows,  the  western  com¬ 
ponent  has  remained  at  a  flux  density  of  approximately 
15%  that  of  the  eastern  component  for  the  past  15  yr. 

Goodman  &  Narayan  (1985)  present  expressions  for  the 
rms  refractive  intensity  fluctuation  as  a  function  of 
density  spectral  index  a  (see  also  Rickett  1986).  In  general, 

~  10%  for  a  <  4  and  >  50%  for  a  >  4.  We  determine 


TABLE  5 

Flux  Density  of  B2  2050  +  36  at  0.61  GHz 


Eastern  Western 

Epoch  (Jy)  (Jy)  Ratio  Reference 


1984  Aug .  3.8  0.6  0.16  Mutel  &  Hodges  1986 

1988  Jun .  3.9  ±  0.1  0.6  ±  0.1  0.15  ±  0.03  FM93 

2000  Jul .  2.864  ±  0.006  0.374  ±  0.014  0.13  ±  0.01  This  work 
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an  observed  value  ^  1.5%  from  the  flux  density  ratios  in 
Table  5,  suggesting  a  <  4.  For  ot  <  4,  the  strong  seattering, 
i.e.,  large  SM,  determined  from  the  observed  scattering 
diameter  (Table  3)  means  that  the  rms  refractive  intensity 
fluctuation  will  be  only  a  slowly  changing  function  of  a. 
From  our  observed  value  of  m^,  we  find  (eq.  [3.1.5]  of 
Goodman  &  Narayan  1985)  a  <  3.85,  consistent  with  the 
value  found  by  fitting  the  angular  diameters  (Table  3). 

3.2.3.  Comparison  with  PSR  B2020  +  28 

The  pulsar  PSR  B2020  +  28  (8?7  from  B2  2050  +  36)  has 
been  the  subject  of  many  monitoring  programs  designed,  at 
least  in  part,  to  measure  refractive  scattering  along  the  line 
of  sight  to  it.  In  all  cases,  it  is  found  that  diffractive  scat¬ 
tering  dominates  refractive  scattering,  as  we  find  for  the  line 
of  sight  to  B2  2050  +  36. 

Smith  &  Wright  (1985)  used  dynamic  spectra  obtained  at 
0.41  MHz  to  measure  the  frequency  drift  rate  and  estimated 
that  the  diffractive  scattering  angle  was  at  least  5  times  that 
of  the  refractive  scattering  angle.  Using  a  structure  function 
analysis  of  a  0.43  GHz  flux  density  monitoring  program, 
LaBrecque,  Rankin,  &  Cordes  (1994)  find  diffractive  scat¬ 
tering  to  dominate  over  refractive  scattering.  Bhat,  Gupta, 
&  Rao  (1999)  monitored  the  pulsar’s  flux  density  at  0.33 
MHz  and  estimated  the  diffractive  and  refractive  scattering 
angles  toward  this  pulsar.  They  found  that  the  diffractive 
scattering  angle  was  roughly  10  times  larger  than  the 
refractive  scattering  angle.  They  also  estimated  a  density 
spectrum  spectral  index  of  a  =  T54.  Stinebring  &  Condon 
(1990)  were  unable  to  observe  any  refractive  modulations  of 
the  pulsar’s  flux  density  in  their  multifrequency  monitoring 
program,  which  is  consistent  with  diffractive  scattering 
dominating  on  this  line  of  sight. 

4.  INTRINSIC  STRUCTURE  AND  KINEMATICS 
OF  B2  2050  +  36 

A  key  assumption  in  our  analysis  of  the  differential  image 
wander  is  that  there  is  little  intrinsic  change  in  the  separa¬ 
tion  of  the  components.  Given  that  on  milliarcsecond 
scales,  many  sources  show  evidence  of  jets  and  bulk  motion, 
this  assumption  clearly  must  be  justified. 

We  have  compared  our  1.67  GHz  observations  with 
those  of  Phillips  &  Mutel  (1981).  This  comparison  is  some¬ 
what  difficult  because,  in  general,  different  antennas  at  dif¬ 
ferent  locations  were  used  in  the  two  VLBI  arrays.  Thus, 
this  comparison  is  necessarily  crude,  but  it  will  be  sufficient 
to  show  that  our  assumption  is  reasonable. 

Phillips  &  Mutel  (1981)  cite  a  separation  of  60  mas 
between  the  two  components  from  observations  in  1979 
September.  If  we  lit  our  1.67  GHz  observations  with  only 
two  components,  we  find  a  separation  of  58  mas,  less  than 
that  of  Phillips  &  Mutel  (1981).  Even  allowing  for  differ¬ 
ences  in  arrays  and  uncertainties  in  the  fitting,  we  conclude 
that  any  intrinsic  change  in  the  separation  between  the  two 
components  is  probably  less  than  1  mas.  Over  the  20.8  yr 
separating  these  measurements,  any  intrinsic  motion  is 
probably  less  than  0.05  mas  yr“^.  At  a  redshift  of  0.354  for 
B2  2050  +  36  (de  Vries  et  al.  2000)  and  for  an  assumed 
Hubble  constant  of  65  km  s“^  Mpc“^,  this  proper  motion 
corresponds  to  a  projected  linear  separation  velocity  of  c. 
Our  limit  is  similar  to  that  found  by  Shaffer,  Kellermann,  & 
Cornwell  (1999)  for  CTD  93,  a  source  with  a  similar 
compact  double  structure. 


5.  CONCLUSIONS 

We  have  conducted  a  multifrequency  suite  of  obser¬ 
vations  on  the  source  B2  2050  +  36.  We  have  used  these 
observations  to  confirm  and  extend  previous  conclusions, 
principally  those  of  FM93,  regarding  interstellar  scattering 
along  this  line  of  sight. 

Our  0.61  MHz  observations  form  the  third  epoch  of 
VLBI  observations  at  this  frequency.  Our  observations  are 
separated  by  12-16  yr  from  the  previous  two  epochs.  As  do 
previous  observations,  we  find  the  structure  of  B2  2050  +  36 
at  this  frequency  to  be  dominated  by  two  components 
approximately  56  mas  apart.  Within  the  uncertainties  of  the 
various  measurements,  this  separation  has  remained 
unchanged  for  the  past  16  yr.  We  can  place  a  reasonably 
strict  limit  of  4  mas  on  any  differential  image  wander  caused 
by  refractive  interstellar  scattering.  In  turn,  this  limit 
implies  that  the  electron  density  power  spectrum  along  this 
line  of  sight  has  a  spectral  index  near  the  Kolmogorov 
value,  with  a  value  of  4  being  highly  unlikely.  Since  the 
diffractive  scattering  diameter  of  the  source  at  this  fre¬ 
quency  is  21  mas,  we  conclude  that  diffractive  scattering 
dominates  along  this  line  of  sight. 

The  most  significant  source  of  uncertainty  in  this  limit  is 
probably  the  velocity  at  which  the  scattering  material  is 
drifting  past  the  line  of  sight.  A  “  high  ”  velocity  of  100  km 
s“^  would  mean  that  the  separations  measured  at  the  three 
epochs  are  independent.  A  “low”  velocity  of  10  km  s“^ 
would  mean  that  none  of  the  epochs  are  independent.  Even 
if  the  velocity  of  the  scattering  material  is  low,  the  separa¬ 
tion  between  our  observations  and  the  previous  ones  is 
probably  long  enough  that,  if  refractive  scattering  were 
important  on  this  line  of  sight,  some  change  in  the  apparent 
separation  of  the  components  would  be  expected.  Nonethe¬ 
less,  future  observations  of  B2  2050  +  36  at  0.61  GHz  or 
similar  observations  of  other  sources  on  lines  of  sight  for 
which  the  refractive  timescale  could  be  estimated  would  be 
useful. 

Our  observations  are  in  good  agreement  with  a  variety  of 
observations  on  the  pulsar  PSR  B2020  +  28  (line  of  sight  8?7 
from  B2  2050  +  36).  These  observations  suggest  that  diffrac¬ 
tive  scattering  also  dominates  refractive  scattering  along  the 
line  of  sight  to  the  pulsar. 

Our  quasi-simultaneous  observations  range  from  0.33  to 
8.4  GHz.  We  have  used  these  multifrequency  measurements 
to  constrain  the  frequency  dependence  of  the  angular  diam¬ 
eter  of  one  component  of  B2  2050  +  36.  At  frequencies 
below  1  GHz,  the  structure  of  the  eastern  component  of  B2 
2050  +  36  is  dominated  by  scattering.  Fitting  for  both  the 
density  spectrum  spectral  index  a  and  the  scattering  diam¬ 
eter  we  find  a  =  3.52  with  the  value  of  4  again  highly 
unlikely.  At  a  nominal  frequency  of  1  GHz,  we  find  the 
scattering  diameter  to  be  0^  =  7  mas.  We  also  place  a  con¬ 
servative  upper  limit  on  the  density  spectrum  inner  scale  of 
li  ^  2000  km. 

We  also  compared  our  1.67  GHz  observations  to  those  of 
Phillips  &  Mutel  (1981),  acquired  approximately  21  yr  pre¬ 
viously.  We  find  no  evidence  for  any  change  in  the  separa¬ 
tion  at  this  frequency.  This  places  a  limit  on  the  projected 
linear  separation  velocity  of  the  two  components  of  c. 

The  National  Radio  Astronomy  Observatory  is  a  facility 
of  the  National  Science  Foundation  operated  under  coo¬ 
perative  agreement  by  Associated  Universities,  Inc.  Basic 
research  in  radio  astronomy  at  the  Naval  Research  Labor¬ 
atory  is  supported  by  the  Office  of  Naval  Research. 
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